Flux Line Channels and Reduced Hall Effect in Impure 2D Superconductors 
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We have studied the dynamics of flux-lines in an impure two-dimensional superconductor using 
numerical Brownian dynamics. In the dirty pinning limit and for small center-of-mass velocities the 
flux lines prefer traveling in narrow channels uniquely determined by the pinning configuration. At 
a certain critical drag force there is a dynamic transition from this glassy state into a more ordered 
crystalline state, in which the channel structure vanishes. The channels are studied by monitoring 
the particle trajectories and the spatial distribution of vortices. By explicitly including a Magnus 
force in the Langevin equations of motion we induce a Hall field and we find that the effective Hall 
angle is strongly reduced by the formation of channels, i.e. the channels introduce a transverse 
energy barrier that guides the motion, predominantly in a direction parallel to the external force. 
The transverse depinning force is larger than the depinning force in the direction of the drag and 
thus the Hall field actually vanishes at finite velocities. We propose that the Hall angle can serve as 
an order parameter for experimental detection of the glass to crystal transition. Our results are in 
excellent agreement with recent predictions made by Giamarchi and Le Doussal. M 
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The flux line lattice (FLL) of thin film and high-T c su- 
perconductors has been found to exhibit many interesting 
properties of which the dynamic phase diagram involving 
both plastic, crystalline, liquid and glassy states is one. 
[Q| In fact this complexity is to be expected since there 
in the system are four strongly competing energy scales; 
thermal fluctuations, vortex-vortex interaction, vortex- 
pinning interaction and the external field. In addition 
the impurities can be realized in many different ways by 
varying the density, distribution, range and strength of 
the interaction potential, in order to model both natural 
and artificially created impurities. 

Much effort has been put into characterizing the nature 
of the dynamic phase diagram and the dynamic phase 
transitions. This is well motivated effort since to 

some extent the energy dissipation of the dynamic vortex 
transitions is resposible for the actual breakdown of the 
superconducting state. 

Koshelev and Vinokour reported computer simulations 
where they monitored the hexagonal order parameter in 
order to characterize the phase diagram spanned by the 
external drag force and the temperature in the presence 
of pinning. || Several attempts have been made to em- 
ploy perturbative methods in order to compute the cen- 
ter of mass velocity as a function of the external force, 
but recently Giamarchi and Le Doussal pointed out that 
perturbation theory in the inverse velocity and disorder 
breaks down for these systems because some modes of 
the disorder are not affected by the motion. In ad- 
dition they also predicted that the vortices should pre- 
fer traveling in elastic channels uniquely determined by 
the disorder configuration. To some extent this predic- 
tion was confirmed already before is was made by the 



experiments of Trauble and Essmann. || Indications of 
channel-like motion has also been seen in computer sim- 
ulations of Brass et. al. but then for a system with an 
interaction potential (Gaussian) far from the true long- 
range potential and at very low pinning densities. || 

In this letter we will use the dynamic phase diagram 
presented by Koshelev et al. [0 as a starting point and 
inspired by the predictions of Giamarchi et al. |l ] we carry 
on the analysis further. 

The starting point is the over-damped (massless) 
Langevin equation for a vortex i, Q,^] 



Vi(f) = F d - VC/ OT (ry) - VU vp {r lk ) 

j^i k 

- avi(t) x z + Ci, 



(1) 



where U vv is the repulsive vortex-vortex potential, U vp 
the attractive vortex-pin potential and Q is a noise term 
satisfying (Ci,a(t')Cj,a'(t)) — 2TSijS a , a >S(t — t') where 
a,a' refer to Cartesian components and T is the tem- 
perature. Frf is an external drag force which corresponds 
to the applied current in the experimental situation and 
a is a coefficient that governs the strength of the Magnus 
force (i.e. the fourth term on the l.h.s. is the Magnus 
force) . 

The equations of motion are simulated using the 
method originally described by Ermak |9) and later used 
for vortices by Brass et. al. |t]], using a numerical time 
step, At = 0.1. 

The vortex system is simulated at intermediate mag- 
netic fields, for which the two-dimensional vortex inter- 
action potential is purely logarithmic, U vv (r) = lnr/ao, 
p"T| where ao is the unit measure of length here taken to 
be the lattice constant of the perfect FLL. This choice of 



1 



potential also defines the unit of energy. In practice we 
use the Ewald summation technique, periodic boundary 
conditions and the minimum image convention in order 
to properly handle the long range forces as described by 
e.g. Frey ct. al. pJ] The interaction potential used here 
is somewhat more ambitious than the potential used by 
Koshelev et al. which was a force shifted logarith- 
mic potential explicitly cut off at 3ao (i.e. not truly 
long-range). Actually a potential with shorter interac- 
tion range more or less just corresponds to a lower mag- 
netic field in the real system. This is due to the pecu- 
liar length scale invariant property of the logarithm func- 
tion. With a full logarithmic interaction all N V (N V — 1)/2 
pair interactions must be taken into account. Thus the 
present simulations are very much more computationally 
demanding than would they be if using a short range po- 
tential and neighbor lists. In practice it however turns 
out that the precise choice of interaction potential makes 
little difference. 

The pinning is realized from a set of N p ran- 
domly distributed attractive Gaussian potentials U vp — 
—A v e~ {j l rv ^'. The parameter set we will consider is: 
A p = 0.006, r p = 0.2, N p = 9500, N v = 20 x 24 = 480, 
0.005, < Fd < 0.25 and for the Magnus coefficient we 
choose the value a — 0.2, i.e. we expect the Hall angle to 
be $h « 11.31° when the drag force is large enough for 
the impurities to be neglected (or for a pure system). 

The schematic phase diagram sketched in Fig. |l| is by 
now well established. [0 
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FIG. 1. The dynamic phase diagram of the vortex system. 
The dashed line starting at T = 0.005 and Fd = is studied 
in the simulations. The vertical dotted line at T = 0.007 
is the absolute melting transition. The lower horizontal 
(dotted) curve displays the plastic-elastic crossover and the 
upper curve (solid) shows the transition from from glass-like 
state into the ordered crystal state. 

By varying the drag force Fd, in our simulation we will 
follow the vertical line starting at Fd — and T = 0.005. 

As also was noted by Koshelev et al. Q the transition 
from the glassy state into the more ordered crystalline 



phase is sharp. The transition can actually be seen di- 
rectly in the vortex pair distribution function which is 
plotted in Fig. |. The set of curves with smaller am- 
plitudes correspond to drag forces below the transition 
whereas the set of curves with large amplitudes corre- 
spond to forces above the transition. Using this some- 
what crude order parameter we determine the transition 
point to be located at Fd « 0.062 for this particular pin- 
ning configuration. 
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FIG. 2. The dynamic glass to crystal transition is clearly 
and intuitively visualized by the vortex-vortex pair distribu- 
tion functions. The set of curves with smaller amplitude cor- 
respond to the drag forces F d = 0.03, 0.04, 0.05, 0.06, below 
the transition, and the curves with larger amplitudes are for 
Fd — 0.07,0.08,0.1 above the transition, respectively. The 
transition occurs at, Fd ~ 0.062. 

The formation of channels is clearly illustrated in 
Figs.|| a-d where every 175th configuration out of totally 
100000 simulation steps is plotted. This gives a time av- 
eraged map of the spatial distribution of vortices in the 
simulation cell. Each dot represents a vortex and thus the 
darkest areas are those most favored by the the flowing 
vortices. The same pinning configuration but different 
starting conditions were used to produce these plots. 

Starting at small forces, F d = 0.002 (Fig.|a), we find 
that most vortices are pinned to the background and the 
total flow is almost completely suppressed by the pinning 
(see also Fig.[|). The few vortices that actually do move, 
creep in channels between the pinned chunks of the lat- 
tice. For Fd — 0.015 (Fig.||b). the channel structure is 
clear. Very few vortices are pinned for long times but the 
mobility is still strongly reduced by the pinning. More 
close to the transition, for Fd — 0.05, the channel struc- 
ture is still present but the vortices have a tendency to 
"tunnel" between the channels. Finally for Fd — 0.08, 
above the transition, only a very weak channel structure 
can be discerned. 
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FIG. 3. Time averaged density maps show the pre- 
ferred spots for a vortex and illustrates how the vor- 
tices prefer traveling in elastic channels. Figures a) - d) 
(from top to bottom) correspond to different drag forces, 
F d = 0.002, 0.015, 0.050, 0.080. 

It is clear that the channel structure is indeed unique 
for each pinning configuration since runs with very dif- 



ferent initial conditions eventually end up with the same 
flow patterns for a given configuration of pins. 

The dependence of the transverse and perpendicular 
center-of-mass velocities with external force (i.e. the IV- 
characteristics) is plotted in Fig. ||. In the same diagram 
also the effective Hall angle is plotted and it is clear that 
the Hall field becomes strongly reduced as soon as the 
channels start to form. Interestingly enough the angle 
is slightly reduced even at velocities far into the "Mov- 
ing crystal" region of Fig.|l|. These velocity versus force 
curves have been averaged over 15 or more realizations 
of the disorder. 
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FIG. 4. The average center of mass vortex velocities in 
the transverse (broken) and horizontal (solid) directions are 
plotted as a function of the drag force Fd . The corresponding 
Hall angle shows that the transverse flow is strongly reduced 
at low velocities due to the presence of the elastic channels. 

The strong dependence of the Hall angle with the drag 
force should be encouraging for designing experiments 
since it may serve as an order parameter for the glass to 
crystal transition. It would also be interesting to make 
more quantitative comparisons with the real system by 
applying the emperical Corter-Casimir formulas to the 
model used here, in order to obtain the correct tempera- 
ture and magnetic field dependencies. [fl3||7l It has been 
argued that the vortex lattice tries to minimize it's power 
dissipation by reorienting in the direction of the drag 
force. Q| This has not yet been verified neither experi- 
mentally nor in simulations although the effect probably 
is of great importance. 

The glass to crystal transition is very sharp in the 
structural and geometrical order parameters and should 
therefore be accompanied by a strong buildup of wash- 
board modes, i.e. harmonic components in the center- 
of-mass velocity. We have preliminary simulational data 
that verify that not only the fundamental mode but also 
higher order frequencies are important in the spectrum of 
the center-of-mass velocity. Jjl| Coupling of these modes 
may enhance the strength of the transition and drive the 
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system into a more ordered state. This and many other 
issues will be dealt with in a forthcoming publication. 




In conclusion we have shown that the flux-lines of a 
2D superconductor at low velocities flow in elastic chan- 
nels that strongly suppresses the Hall field and that the 
Hall angle vanishes already at finite velocities. The flow 
patterns are uniquely determined by the pinning config- 
uration. 
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